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^) Abstract. We derive a simple relationship between the energy emitted in 

gravitational waves for a narrowband source and the distance to which that emission 
can be detected by a single detector. We consider linearly polarized, elliptically 
polarized, and unpolarized gravitational waves, and emission patterns appropriate for 
each of these cases. We ignore cosmological effects. 
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pg The sensitivity of a gravitational- wave detector to a transient signal ("burst") is 

O reasonably well characterised by the expectation value of the matched-filter signal-to- 

^ noise ratio (SNR) of the burst \T\. This expectation value can be computed as an integral 

CO of the burst signal spectrum divided by the detector noise spectral density. For the most 

commonly considered gravitational- wave transient signals, the inspiral of neutron star 
binaries, this leads to the well-known average (or "sensemon"E]) range [2] giving the 
KS effective volume within which the expected SNR in a single detector would be above 

some threshold pdet'- 

Here A^ is the chirp mass of the binary and f^/s = J dff "^^'^S ^ where 5* is the noise 
power spectrum of the detector. The factor 1.77 arises from integration over possible 
positions and orientations of the binary. The average range T^bns is defined such that 
for a homogenous isotropic distribution of sources with rate density Af the mean rate of 
detections will then be 

iV=^7^(7^BNs)'A^. (2) 

J The other commonly used range for inspirals is the "horizon" range, defined as the maximum distance 
at which an optimally positioned and oriented binary would produce an expected SNR of at least pdct ■ 
The horizon range is a factor 2.26 larger than the average range. 
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The average range is clearly a useful tool for assessing the scientific capabilities 
of a detector. Unfortunately, while the expected gravitational-wave emission of binary 
systems is well-understood, gravitational wave transients from most other potential 
sources such as supernovae, long gamma-ray bursts, and soft gamma repeaters are not 
known reliably. In this technical note we present a simple definition of the effective 
sensitive range of a gravitational-wave detector to generic bursts, based on the total 
energy emitted in the burst and the peak frequency of the emission. This range is 
appropriate for a variety of possible signal polarizations and source emission patterns, 
and requires only that the signal bandwidth is smaller than the frequency range over 
which the detector noise spectrum varies significantly. 

We begin by determining how the total energy -Egw carried in a gravitational-wave 
burst is related to the measure of signal strength commonly used in burst searches, the 
root-sum-square amplitude h^ss- We then relate these measures to the expected signal- 
to-noise ratio p, and define an average range analogous to equation ([I]). Finally, we 
compare this range to the results of recent LIGO- Virgo searches for gravitational-wave 
bursts [31 HE]. 

2. Relating -Eqw to hj-gg 
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We first relate the total energy emitted in gravitational waves, -Eqw; to the LIGO- Virgo 
standard measure for burst amplitude at the detector, hj-ss- The latter is defined by 

hr..= I dt [hl{t) + hlit)] (3) 

df \\KU)? + \hAf)?\ ■ (4) 

The flux (energy per unit area per unit time) of a gravitational wave is 

Fgv, = Tiy^{hl{t) + hl{t)) , (5) 

iovrG 

where the angle brackets denote an average over several periods. For a burst of duration 

< T we can compute the average by integrating over the duration: 

^3 1 /•T'/2 



-^GW 



IQttGT 
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-T/2 
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^ , dt 
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+ (same. 



hl{t) + hl{t)_ 

df^'-f''{i27^f')hX{f') 



-)■ X 



(6) 

dfe-'^-f\-i2Txf)h^{f) 



(7) 



Since /i+^x -^ outside —T/2 < t < T/2, we may extend the time integration to 
t — 7- ±oo. The time integral then evaluates to a delta function, 5{f — /'), giving 

_ TIC" 1 

4GT 



GW 



dffi\K{f)\^ + \K{f)\- 



(8) 
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2.1. Isotropic emission 

To compute the total energy -Egw emitted, we need to integrate the flux Fqw assuming 
some emission pattern. Let us first assume isotropic emission. Then for a source at a 
distance r 



EGw = 47rr2TFGw (9) 

r' rdff{\hM)\' + \hM)\') ■ (10) 



TT^C^ 



G 

If we assume that the signal is narrowband with central frequency /o, we obtain 






^GW = ^r%Xss- (11) 

2.2. Linear motion emission 

Axisymmetric motion will produce linearly polarized emission with pattern 

h+{t) = sm^{L)h{t), (12) 

/ix(t) = 0, (13) 

where l is the angle between the symmetry axis and the line-of-sight to the observer, 
and we have selected a polarization basis aligned with this symmetry axis. The energy 
emitted in a narrowband signal is then 

EGw = j^r' d{cosi) dX / dff{sm\i)\hif)\'- 

o — - _L v U </ — OO 



2^3 



TT C 



^^ ^ rVo'/^L, (14) 

where A is the azimuthal angle in the source frame. This is 8/15 times the result for 



isotropic emission, (11). 



Note that in writing (14) we have defined /ij-ss as the root-sum-square amplitude 
from an optimally oriented source {l = tt/2 in this case). This differs slightly from 
the standard LIGO- Virgo definition, which includes the inclination factors. In practice, 
however, LIGO-Virgo papers to date have typically simulated optimally oriented sources. 

2.3. Rotating system emission 

Rotational motion (such as from a circular binary) will produce emission with pattern 

h+{t) = -{1 + cos'^it)) A{t) cos <l>{t) , (15) 

/ix (t) = cos(0 A{t) sin ^{t) , (16) 

where l is the angle between the rotation axis and the line-of-sight to the observer, 
and we have again selected a polarization basis aligned with this symmetry axis. We 
assume A{t) varies slowly enough compared to $(t) that h^^ and hx are approximately 
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orthogonal. This produces an eUiptically polarized signal at the detector. The energy 
emitted in a narrowband signal is 

= ~r^fXss, (17) 

where h{f) is the Fourier transform of A{t) cos<l>(t) and we have again used /irss for an 
optimally oriented source (t = 0). The expression for energy emitted is 2/5 times the 



result for isotropic emission, (11). 



3. Relating -Eqw to Signal- To-Noise Ratio 

The detectability of a generic signal is determined mainly by its expected signal-to- 
noise ratio p for a matched filter. (The time-frequency volume Vtf of the signal is also 
important when Vtf ^ 1 ID E|)- For a narrowband signal, p has a simple relationship 
to the hj-ss amplitude. We start from 

where S{f) is the one-sided noise power spectrum, and F+_x(6',0, ■?/') are the antenna 
responses to the sky position [6, (p) and polarization ip of the gravitational wave. We 
may expand the square in (18) and drop the h^h*^ terms for most signals of interest: 
for eUiptically polarized signals the two waveforms are orthogonal, while for linearly 
polarized signals /ix = 0. The waveforms are also orthogonal in the unpolarized case, 
where the two polarizations are independent stochastic timeseries. An example is white- 
noise bursts [5j. Assuming a narrowband signal, we find 

P^ = e^^, (19) 

where we define the angle factor 

{Fl{9, 0, ijj) + F^ {9, 0, ip) isotropic unpolarized 

F2(0,0,^)(l±Hf^)2 + F2(0,0,7/')cos2(O elliptical (20) 

Fl{9,(f),ip)2 sin"^ l linear 

Note that all dependence on the four angles 9, 0, ip, and l is contained in 9. Substituting 



(11), (14), or (17) gives 



where a = 1 for isotropic emission, 8/15 for linearly polarized emission, and 2/5 for 
circularly polarized emission. 
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4. Effective Range 

We can now combine the results for -Egw and p to compute tlie typical distance to which 
a source is detectable. We will follow the approach used in Section V of |2]. 

Consider a homogenous isotropic distribution of sources with rate density A/". A 
signal from a given source will be detectable if the received signal-to-noise is above some 
threshold value pdet- The mean rate of detections will then be 



N = AttX / drr^P{p'^ > pi 



ct^ 



(22) 



Here P{p^ > p^et) i^ the probability that the signal-to-noise of a source at given distance 
r with random 6, (^, ip, and t will be above threshold. Using (21), we may write this 
probability as 

Pip' > pL) = PiQ' > r'/rl) , (23) 

where we have defined the fiducial distance 
2 G -Egw 



° avr2c3 S{fo)f^pl 
Our detection rate is thus 

4 



et 



N, 



det 



-Tcr^N 



3 / dxx'P{Q' >x' 



The integral is easily evaluated numerically: 

0.0978 unpolarized 

dxx^P{Q'^ > a;^) = <j 0.0287 elliptical 

0.0537 linear 



(24) 



(25) 



(26) 



Following [2], we define the effective detection range TZeS as the radius enclosing a 
spherical volume V such that the rate of detections is AfV: 
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(27) 




(28) 


0.665 unpolarized 




0.698 elliptical 


(29) 


0.745 linear 





We note that for all three cases (unpolarized, linear, and elliptical polarizations), /3 is 
equal to l/-\/2 to within a few percent. A convenient approximation is thus 



7^ 
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E, 



GW 



27r2c3 5(/o)/oVi 



ct 



1/2 



(30) 



With this definition the mean rate of detections for a homogenous isotropic distribution 
of standard- candle (fixed -Egw? /o) burst sources with rate density Af is 

iV=^7r(Kff)'Ar. (31) 
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5. Example: LIGO-Virgo Science Runs, 2005-2010 



As an example, we apply our effective range formula (30) to the LIGO- Virgo network 
during their 2005-07 and 2009-10 science runs. The results of the search for generic 
gravitational- wave bursts are reported in |31 SI E]- Approximately 1.8 yr of coincident 
data were analysed from the three LIGO detectors (HI, H2, LI) and the Virgo detector 
(VI). No gravitational waves were detected, and limits were placed on the rate, 
amplitude, and energy content of gravitational waves. 

Figure [l](a) shows an example noise spectrum from each of the detectors that 
participated in the 2009-10 run (data obtained from [7j). The HI detector had the lowest 
noise level across most of the search frequency band, so for convenience we use its noise 
spectrum S{f) for our range calculations. The other quantity required for defining the 
range is the SNR threshold pdet, which is the threshold at which the detection efficiency 
is 50%. Comparing the h^ss amplitude limits for linearly polarized sine-Gaussian bursts 
and unpolarized white-noise bursts (Tables II and IV of [5J, Table II of [4J, Fig. 3 of 
[3]) show that they correspond to pdet — 20 to 30 as measured against the HI S6 noise 
spectrunjH depending on the waveform, network, and data set. Since we expect the rate 
limits to be dominated by the most sensitive data (due to volume scaling), we select 
Pdet = 20 for our estimates. 

Figure mb) shows the effective range (30) predicted assuming pdet = 20 and the HI 
noise curve smoothed to 10 Hz resolution. The left-hand scale (Mpc) assumes -Egw = 
10^^ MqC^, which is the approximate maximum gravitational- wave emission possible 
from long gamma-ray bursts under the most optimistic scenarios [HI [9], [TOl [HI [121 [IS [E] ■ 
The right-hand scale (kpc) assumes -Egw = 10~^MqC^, which is a typical energy emission 
in simulations of core-collapse supernovae fibl [16] . The maximum ranges in the two 
cases are approximately 10 Mpc (10 kpc) for signal frequencies around 100 Hz - 200 Hz, 
dropping to below 1 Mpc (Ikpc) by 1000 Hz. 

Figure [l](c) shows the -Egw predicted by (30) to be required for a source at a fixed 



distance of 10 kpc to produce an expected SNR equal to pdet = 20. The dots are the 
actual -Eqw values for a variety of waveforms and the 2009-10 HILIVI network, as 
reported in Fig. 7 of [S|. Figure [l](d) shows the 90% confidence rate density limit (rate 



per unit volume) predicted by (31 ) for a homogeneous isotropic distribution of standard- 
candle sources with Eqw = IMqC^, assuming no detections in the 2005-07 and 2009-10 
searches (so that A^ < 2.3 at 90% confidence). The dots are the approximate rate 
density limits for linearly polarized sine-Gaussian waveforms set by the combined 2005- 
07 and 2009-10 data sets (Fig. 6 of [SJ). Despite the fact that we use a single sample 
noise spectrum and SNR threshold to represent all networks and both science runs, the 
predicted limits are a reasonably good match to the measured limits in each case. 



§ The SNR threshold p^Gt can be estimated from the amphtude Hmit h^.g^" using (19 20 1: pdct 
0rms^rs°^Vv^'5'(/o) where Ornis = ((^+ + ^x )e,0,i/.)^^^ = ^2/5 for optimally oriented sources. 
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Figure 1. (a) Example noise spectrum from eacli of tlie detectors in the 2009-10 
LIGO- Virgo science run. (b) Effective range predicted by (30) using pdct — 20 and the 



HI noise curve smoothed to 10 Hz resolution. The left-hand (right-hand) scales assume 
^GW = 1O~^M0C^ (1O~^M0C^). (c) i?GW predicted by (30) to be required for a source 
at a fixed distance of lOkpc to produce an expected SNR equal to pdot = 20. 



The 



dots are the actual -Egw values for a variety of waveforms and the 2009-10 HILIVI 



network [5]. (d) 90% confidence rate density limits predicted by (31 1 for a homogeneous 
isotropic distribution of standard-candle sources with E'gw — IMqc^. The dots are 
rate density limits for linearly polarized sine-Gaussian waveforms set by the combined 
2005-07 and 2009-10 data sets 0. 
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